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forward  scattering  by  elastic  spherical  shells  and  rigid/soft  spheres,  for  nondimensional  frequency  values 
20  <  ka  <  80,  and  scattering  angles  0°  <  0  <  10°.  There  arc  only  small  differences  among  the  Fraunhofer,  Fresnel, 
and  rigid/soft  predictions  for  differential  scattering  cross  sections.  They  arise  mainly  from  the  slightly  varying 
degrees  of  compression  (in  angular  space)  of  one  pattern  relative  to  another.  The  magnitudes  of  the  maximums 
predicted  by  the  four  methods  are  in  good  agreement,  but  the  minimums  do  not  agree.  Calculations  for  1%,  5%, 
and  10%  spherical  steel  shells  indicate  that  variation  of  shell  thickness  has  only  a  small  effect  on  forward  scattering. 


1.  INTRODUCTION 

Forward  scattering  has  been  used  in  particle  sizing  in  optical*  and  acoustic^  measurements.  Both  Fraunhofer 
diffraction  and  Fresnel  diffraction^  have  been  used  as  predictive  tools  in  calculations  of  scattering  from  objects  of 
suitably  small  sizes.  Since  acoustic  waves  do  penetrate  elastic  objects,  but  the  Fraunhofer  and  Fresnel  diffractions 
do  not  include  elasticity,  the  validity  of  these  physical  optics  approaches  even  in  forward  scattering  may  be  ques¬ 
tioned.  In  the  present  report  we  compare  predictions  of  forward  scattering  from  elastic  spherical  shells  and  rigid/ 
soft  spheres  using  known  exact  formulas,^  to  Fraunhofer  and  Fresnel  predictions. 

The  methods  of  physical  optics  have  been  applied  to  acou.>(ic  scattering  problems  by  many  auihors^’^"'^  with 
varying  degrees  of  success.  The  most  well  known  of  these  techniques  is  the  Kirchhoff  approximation  (Ref.  5  and 
references  therein).  It  assumes  that  the  incident  field  and  its  normal  derivative  are  unperturbed  in  an  open  aperture, 
and  zero  on  the  shadow  side  of  a  scatterer.  Further,  the  outgoing  signal  on  the  insoniHed  side  of  a  scatterer  is  taken 
to  be  the  incident  signal  multiplied  by  its  plane  wave  reflection  coefficient  at  every  point  on  the  scattering  surface. 
The  application  of  these  assumptions  to  backscattering  from  semi-rnfiniie  plates,*  and  from  spheres  and 
cylinders^*  have  shown  that  the  Kirchhoff  method  is  unreliable  in  these  instances. 

Applications  of  the  Kirchhoff  method  to  forward  scattering  problems  have  been  more  successful.  Forward  scat¬ 
tering  from  a  semi-infinite  plate  has  been  discussed  by  Officer.’  Pierce  and  Hadden  find'”’"  that  for  small  forward 
angles,  the  Kirchhoff  method  provides  a  good  approximation  to  scattering  by  a  wedge.  This  result  is  independent 
of  the  wedge  angle.  Hunter,  Lee,  and  Waag^  have  experimentally  measured  forward  scattering  patterns  of  single 
nylon  filaments.  They  demonstrated  that  the  patterns  were  well  explained  by  a  Huygens  construction,  calculated 
numerically.  The  Fraunhofer  and  Fresnel  calculations*  discussed  below  arc  subsets  (subsets  that  contain  all  the 
important  elements)  of  the  Huygens  construction.  In  an  earlier  work,'*  this  author  compared  the  total  field  due  to 
an  infinite  rigid  cylinder  at  small  forward  angles  calculated  using  exact  formulas,'*  to  that  calculated  using  the 
Kirchhoff  method  adapted  from  Officer’  for  an  opaque  strip.  The  two  predictions  were  in  good  agreement.  For 
forward  scattering,  the  Kirchhoff  method  is  identical  to  Fresnel  diffraction.*-’  Since  the  concepts  of  Fresnel  diffraction 
and  Fraunhofer  diffraction  arc  well  established  as  special  cases  of  a  unified  theory  based  on  the  Fresnel-Kirchhoff 
integral,*  we  prefer  to  retain  the  Fresnel  label  instead  of  the  Kirchhoff  label.  We  show  that  in  the  examples 
considered  here,  both  Fraunhofer  diffraction  and  Fresnel  diffraction  provide  good  approximations  to  the  scattered 
field  as  well  as  the  total  field  due  to  clastic  spherical  shells  and  rigid/soft  spheres.'^ 
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2.  EXACT  FORMULAS 


The  well-known  exact  formulas*  for  spherical  elastic  and  rigid/soft  scatterers,  are  summarized  in  the  present 
section.  For  spherical  scatterers,  the  complex  pressure  of  the  incident  plane  wave  can  be  written  as 

Pine  =  exp  (ikreosG) , 

where  k  is  the  wave  number,  r  is  the  distance  of  the  observation  point  from  the  center  of  the  sphere,  and  6  is 
spherical  polar  angle,  with  6  =  0  for  forward  direction.  Similarly  for  the  scattered  signal, 


P«.u  =  I  (2n  +  OA^PJeosG)  =  f-p]exp(ikr)f(n).  (2) 

Here,  (cos  G)  are  the  Legendre  polynomials,  and  the  coefficients  A„  are  determined  from  elastic  boundary 
conditions.*  The  formulas  are  for  the  asymptotic  case  kr  »  1.  For  rigid/soft  boundary  conditions,* 

A„  (rigid)  =  -  j'„(ka)/h‘„(ka),  and  A„  (soft)  =  -  j„(ka)/h„(ka) .  (3) 

Here  the  j„  are  the  spherical  Bessel  functions,  and  the  h„  are  the  spherical  Hankel  functions  of  the  first  kind. 


The  differential  scattering  cross  section  in  any  direction  Q  is  given  by 


(4) 


The  total  field  at  any  point  is  given  by  the  coherent  sum  of  incident  and  scattered  fields  as 

Pioi  “  Pine  Pscaii  ■ 


3.  PHYSICAL  OPTICS  FORMULAS 


In  the  physical  optics  approach  to  the  forward  scattering  problem,  we  replace  the  scattering  object  (sphere)  with 
an  opaque  screen  having  the  same  area  of  cross  section.  Consider  the  case  when  the  wavefront  due  to  the  incident 
signal  coincides  with  the  screen  (Fig.  1).  We  assume  the  signal  and  its  derivative  on  the  shadow  side  of  the  screen 
to  be  zero,  as  in  the  Kirchhoff  method.  The  signal  at  every  other  point  on  the  wavefront  is  assumed  to  be  unperturbed. 
The  complex  pressure  at  any  forward  point  is  then  given  by  the  Fresndl-Kirchhoff  integral^ 

p,o.  =  C  n  exp  (ikf  (^,  Ti)]d^dTi  (6) 

insonifted 

where  k  is  the  wave  number,  the  function  f(^,Ti)  =  ((r-r')  +  (s-s')]  (Fig.  1),  and  the  three-dimensional  Cartesian 
coordinates  of  P',  Q,  and  P  are  given  by  P'(x',y',z'),  Q(4,'n,0),  and  P(x,y,z).  For  incident  plane  waves  of  unit  strength 
and  small  forward  scattering  angles,  the  constant  C  is  given  by^ 


C'  =  -  i 


.  k 


2n 


exp(iks ') 


(6a) 


When  the  coordinates  x',  y’,  ti,  x,  and  y  are  small  compared  to  both  r’  and  s',  we  can  expand  f(^,Ti)  in  a  power 
series  as 


f(^.Tl)  =  - 


x'4  +  y’Ti  x^-^yTl  ^  +  TiVl  1 


+ - . 


(7) 
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Fig.  1.  Source,  receiver,  coordinate  axes,  and  the  position  of  the  wavefront  as  it 
coincides  with  the  screen. 


The  integral  in  Eq.  (6)  runs  over  the  entire  insonified  region,  viz.  the  wavefront  extending  to  +  /  -  “»  in  ^  and  t), 
but  excluding  the  opaque  screen.  (In  the  following  we  only  evaluate  integrals  over  a  finite  region  of  the  wavefront, 
consistent  with  the  condition  of  smallness  of  ^  and  tj  relative  to  r’  and  s'.) 


For  better  physical  insight,  let  us  rewrite  Eq.  (6)  as 

p  =C'  J  Jexp[ikf(^,Tl)]d^dn-C’  |Jexplikf(tTl)ld4*l  =  Pi^ (8) 

-  iite 

where  the  integral  labeled  “disc”  runs  over  the  opaque  screen.  Equation  (8)  can  be  understood  as  a  result  of 
Babinet’s  principle.^  The  first  term  Pi„j  in  Eq.  (8)  represents  the  field  due  to  the  entire  infinite  wavefront,  the  field 
in  the  absence  of  any  scatterer.  The  second  term  p^;^  in  Eq.  (8)  represents  the  field  due  to  the  complementary 
geometry,  an  insonified  opening  in  an  infinite  opaque  screen. 


The  terms  p,o,  and  Pi„e  in  Eq.  (8)  have  exactly  the  same  meaning  as  the  quantities  p,o,  and  Pi„j  in  Eq.  (5)  Therefore, 
the  quantity  (-Pdi,c)  in  Eq.  (8)  corresponds  to  the  term  p„,„  in  Eq.  (5).  If  both  the  diffraction  theory  and  the 
scattering  theory  make  identical  predictions,  then  for  the  same  incident  field  Pu,*  in  Eqs.  (5)  and  (8),  (-Pdi,*)  of 
Eq.  (8)  will  be  equal  to  p„,„  of  Eq.  (5)  Following  the  scattering  theory  definition  of  differential  scattering  cross 
section  in  Eq.  (4),  we  define  the  corresponding  quantity  for  diffraction  theory  as 


(9) 
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When  the  quadratic  and  higher  terms  in  f(^,Tl)  in  Eqs.  (6)  and  (7)  can  be  ignored,  we  have  the  case  of  Fraunhofer 
diffraction.^  Bom  and  Woli^  demonstrate  the  condition  for  this  as 


r.s» - - - 

where  X  is  the  wavelength.  When  the  quadratic  terms  are  important,  we  have  Fresnel  diffraction.^ 


(10) 


Let  us  now  discuss  the  calculation  of  Pj^  and  p^  of  Eq.  (8)  under  the  two  approximations.  Evaluation  of  p^ 
using  Fraunhofer  approximation  poses  no  conceptual  difficulty.  Thus  within  the  Fraunhofer  approximation,  using 
some  change  of  variables  and  suaightforward  integration  we  obtain^ 


p  ^  (Fraunhofer)  =  C  ]]  exp 


diic 


[- 


n  ^  x4  +  yq 


d^dq 


2Jj(ka  sin  6) 
ka  sin  6 


(11) 


Here  9  is  the  angular  separation  of  a  point  on  the  diffraction  pattern  measured  from  the  forward  direction,  and  is 
the  same  as  the  angle  9  in  the  previous  section. 

To  calculate  p^ue  using  the  Fresnel  approximation.  Bom  and  Wolf^  choose  a  coordinate  system  such  that  the 
linear  terms  in  f(^,q)  are  identically  zero.  Then 

f(t  11)  =  y(-Jr  +  cos2  8  +  q*)  + - ,  (12) 

where  8  is  the  (small)  angle  between  the  line  FOP  and  the  normal  to  the  opaque  screen.  Defining  new  variables  u 


and  V,  a  new  constant  b’,  and  indefinite  integrals  and  S„v  as 

(13a) 

(13b) 

b'  -  ** 

k(l  /  r'  +  1  /  s’)  cos  8 

(13c) 

Cu,  3  j]cos^y(u^+  v^)jdudv 

(13d) 

S  a //sin 

(13e) 

we  can  write 

PdiK  )  =  C’  b’  //  exppy( u^  +  v7)jdu dv 

disc 

=  -  Y  cxp(iks  •)  [C„,(disc)  +  iS  „,(disc)]  • 

(14) 
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Here,  (disc)  on  the  right  hand  side  of  Eq.  (14)  means  that  the  integrals  run  over  the  opaque  screen.  Using  trigono¬ 
metric  expansions  of  cosines  and  sines  we  obtain 

C^^(disc)  =  2  J  cos(|-u*]qv^_^(u)]du -2  /  sin[yu*]s[v_^Ju)]du  (15a) 

S„,(disc)  =  2  i  sin[|u*)c[v^^(u)]du-H2  j  cos[|u'js[v^_^(u)]du  (15b) 

where  the  line  integral  in  the  variable  u  runs  from  u„u,  to  u„„  (Fig.  2).  The  quantity  v,„„  is  given  by 

r  -|i  n 

''m„  =  La’^-(u  +  u^)  j  (16) 


“■'[^7  4]]  “  <■’> 

is  the  reduced  radius  of  the  sphere  in  u-v  space.  The  quantity  Uq  is  the  (positive)  distance  in  u-v  plane  between  the 
origin  of  the  coordinate  system  and  the  center  of  the  opaque  disc  (Fig.  2).  This  distance  is  also  equal  to  the  distance 


Fig.  2.  The  circular  opaque  screen  (projection  of  sphere) 
in  the  u-v  plane,  as  seen  by  an  observer  moving  in  the 
direction  of  the  incident  signal. 
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in  the  obsen^ation  plane  (also  measured  in  reduced  space)  of  a  point  on  the  diffraction  pattern  from  the  center  of 
the  pattern.  Here  we  have  approximated  cos  d  by  unity,  using  the  fact  that  5  is  very  small.  The  quantities  C  and 
S  are  the  well  known  Fresnel  integrals^  given  by 

C(x)  = /cos^-|u^jdu  (Iga) 

S(x)=  Jsin^yu*^du  .  (18b) 

Formulas'^  to  calculate  C(x)  and  S(x)  are  given  in  Appendix  A.  We  calculated  (disc)  and  (disc)  of  Eqs.  (I4) 
and  (IS)  by  numerical  integration  using  Simpson's  rule. 

The  incident  field  Pi„c  is  a  plane  wave  of  unit  strength,  given  by 

Pine  (Fraunhofer)  =  exp  (iks’  cos  0)  (19) 

in  polar  coordinates  (s',  0),  used  in  the  Fraunhofer  calculation.  As  mentioned  above,  the  Fresnel  calculation  uses  a 
special  coordinate  system^  in  which  the  diffracted  field  is  calculated  on  a  plane  perpendicular  to  the  incident  signal 
direction,  at  distance  s'  from  the  scatterer.  For  this,  the  incident  signal  is  given  by 

Pine  (Fresnel)  =  exp  (iks')  (20) 


4.  NUMERICAL  RESULTS 

Here  we  present  the  results  of  our  calculation  of  differential  scattering  cross  section  and  the  total  field  based  on 
the  Fraunhofer  and  the  Fresnel  approximations.  These  results  are  compared  with  conesponding  results  obtained 
using  exact  formulas^  for  elastic  spherical  shells  and  rigid/soft  spheres.  In  all  examples,  we  have  calculated  the 
diffraction  pattern  in  the  far  field,  with  (s'A)  s  1000. 

Figures  3-S  show  the  differential  scattering  cross  section  da/dQ  plotted  vs.  the  forward  scattering  angle  6,  for 
ka  =  20,  40,  and  80.  There  are  only  small  differences  among  the  Fraunhofer,  Fresnel,  and  rigid/soft  predictions  for 
each  case.  They  arise  mainly  from  the  slightly  varying  degrees  of  compression  (in  angular  space)  of  one  pattern 
relative  to  another.  The  magnitudes  of  the  maximums  predicted  by  the  four  methods  are  in  good  agreement,  but  the 
minimums  do  not  agree. 

Figures  6-8  show  the  squared  magnitude  of  the  total  field  p,M  plotted  vs.  the  forward  scattering  angle  0,  for 
ka  =  20,  40,  and  80.  The  total  field  p,o,  was  calculated  for  an  incident  field  p^  of  unit  strength.  In  each  case, 
agreement  among  the  Fraunhofer,  Fresnel,  and  rigid/soft  predictions  is  remarkably  good.  The  diffraction  pattern 
shows  regularly  spaced  maximums  and  minimums.  In  ail  cases  the  asymptotic  value  (as  6  increases)  is  0  dB, 
consistent  with  the  normalization  of  the  incident  signal  strength  to  unity,  and  the  fact  that  the  scattered  field  strength 
falls  off  rapidly  for  increasing  6.  The  unit  normalization  of  the  incident  signal  also  explains  why  the  oscillations  are 
centered  at  0  dB. 

The  effect  of  shell  thickness  on  forward  scattering  is  illustrated  in  Fig.  9  for  differential  cross  section  do/df),  and 
in  Fig.  10  for  the  total  field  p,^,,  for  spherical  steel  shells  of  thicknesses  1%,  5%,  and  10%.  The  ka  value  is  20 
respectively.  Fresnel  and  rigid/soft  predictions  are  shown  for  comparison.  Shell  thickness  has  only  a  minimal  effect 
in  these  examples. 
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Fig.  6.  Variation  of  the  total  acoustic  field  strength  with  forward 
scattering  angle  6  for  spherical  rigid/soft  spheres,  compared  with 
Fraunhofer  and  Fresnel  diffractions,  ka  =  20. 
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Fig.  7.  Similar  to  Fig.  6,  for  ka  =  40. 
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Fig.  9.  Variation  of  differential  scattering  cross  section  da/df2  with  forward  scattering 
angle  8  for  sphere.  Results  of  1%  ,  S%,  and  10%  steel  shells  are  compared  with  rigid  and 
soft  body  results,  and  Fresnel  diffraction,  ka  =  20. 
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Fig.  10.  Variation  of  the  total  acoustic  field  p^,  with  forward  scattering  angle  9  for 
sphere.  Results  of  1%  ,  5%,  and  10%  steel  shells  are  compared  with  rigid  and  soft  body 
results,  and  Fresnel  diffraction,  ka  =  20. 
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5.  CONCLUSIONS 


Our  numerical  results  demonstrate  that  both  Fraunhofer  and  Fresnel  diffractions  provide  good  approximations 
forward  scattering  by  elastic  spherical  shells  and  rigid/soft  spheres,  for  nondimensional  frequency  value 
20  <  ka  <  80,  and  scattering  angles  0*  <  0  <  10®.  There  are  only  small  differences  among  the  Fraunhofer,  Fresnel 
and  rigid/soft  predictions  for  differential  scattering  cross  sections.  They  arise  mainly  from  the  slightly  varyin' 
degrees  of  compression  (in  angular  space)  of  one  pattern  relative  to  another.  The  magnitudes  of  the  maximums 
predicted  by  the  four  methods  are  in  good  agreement,  but  the  minimums  do  not  agree.  Calculations  for  1%,  5% 

10%  spherical  steel  shells  indicate  that  variation  of  shell  thickness  has  only  a  small  effect  on  forward  scattering 
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APPENDIX  A.  FORMULAS  TO  CALCULATE  FRESNEL  INTEGRALS 
For  X  <  2,55  we  have  used’^ 


n  a  0 


(2n)!(4n  +  1) 


For  x  >  2.55  we  have  used'’ 


„7o  (2n+  l)!(4n  +  3) 


C(x)  =  ^  +f(x)sin^-|x^y-g(x)cos^-|x^^ 
S(x)  =  ^  -f(x)cos(^-yx^j-g(x)sin 


,,  ,  ,  ;  (-  l)"’l.  3....(4m-  1) 

7txf(x)=l+  Z  - 

(tcx^) 


2m 


1  -  (-l)"’l.3....(4m+ 1) 

nxg(x)=  - r  +  z  - 

7tx'‘  ni  =  l 


(ttx^) 


2m  +  1 
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